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We solve for the electronic Raman scattering response functions on an infinite-dimensional hyper- 
cubic lattice employing dynamical mean field theory. This contribution extends previous work on 
the nonresonant response to include the mixed and resonant contributions. We focus our attention 
on the spinless Falicov-Kimball model, where the problem can be solved exactly, and the system can 
be tuned to go through a Mott-Hubbard-like metal-insulator transition. Resonant effects vary in 
different scattering geometries, corresponding to the symmetries of the charge excitations scattered 
by the light. We do find that the Raman response is large near the double resonance, where the 
transfered frequency is close to the incident photon frequency. We also find a joint resonance of 
both the charge-transfer peak and the low-energy peak when the incident photon frequency is on the 
order of the interaction strength. In general, the resonance effects can create order of magnitude (or 
more) enhancements of features in the nonresonant response, especially when the incident photon 
frequency is somewhat larger than the frequency of the nonresonant feature. Finally, we find that 
the resonant effects also exhibit isosbestic behavior, even in the Aig and -B2g sectors, and it is most 
prominent when the incident photon frequency is on the order of the interaction energy. 

PACS numbers: 78.30.-j,71.10.-w,71.27.+a,71.30.+h,78.20.Bh 



I. INTRODUCTION 

Electronic Raman scattering has long been used as a 
direct probe of the charge excitations of different mate- 
rials. Experiments have shown a number of interesting 
phenomena, especially in correlated materials. A mate- 
rial independence for Raman scattering has been seen in 
a number of different correlation gap (insulating) materi- 
als ranging from FeSi^, to SmBg^, to Ca3Ru207'^ to high 
temperature superconductors'*^^'^ The Raman response 
shows a gap opening at low temperature, but with the 
gap about ten times larger than the onset temperature 
where the gap starts to form. In addition, an isosbestic 
point is often seen, where the Raman response at one 
value of frequency is independent of temperature (at low 
temperature), and curves for different temperatures ap- 
pear to cross at a single point. Resonant effects are even 
more interesting, as it is believed that the resonance can 
cause an enhancement of the nonresonant signal by or- 
ders of magnitude, and allow small signals to become ob- 
servable. What remains unknown is whether these reso- 
nant enhancements dramatically change the shape of the 
underlying nonresonant response. 

The theoretical description of electronic Raman scat- 
tering has lagged behind experiment. Recently, dynami- 
cal mean field theory (DMFT) has been employed to cal- 
culate the nonresonant response in the Falicov-Kimball^ 
and Hubbard models* and to examine inelastic x-ray 



scattering as well^. It was found that the theoretical cal- 
culations of the nonresonant response show much of the 
behavior seen in experiment, including the large gap rela- 
tive to the onset temperature and the generic appearance 
of an easily observed isosbestic point in the Big channel. 

However, it is well known that many of the Raman 
signals in correlated metals and insulators display com- 
plicated dependences on the incoming photon frequency 
uji. The resonant behavior of the Big two-magnon fea- 
ture at roughly 340 meV has been well studied in the 
parent insulating cuprates La2Cu04, YBa2Cu306, and 
Sr2Cu02Cl2'', where a resonance is found for incident 
photon energies near 3 eV. Although recent progress 
has been made^°'^*, the reason for this resonance is not 
clear since the resonance frequency lies above the opti- 
cal absorption edge frequency measured in the dielectric 
response^, and the photon energy is much larger than the 
location of the resonance peak in the response function. 

The general question of how the low-energy features 
(such as particle-hole excitations near the Fermi level) 
and high-energy (such as charge-transfer excitations) 
change under resonant conditions remains relatively un- 
explored. Most treatments for Raman scattering in in- 
sulators have focused on only the spin degrees of free- 
dom (Heisenberg limit) in two dimensions. For the case 
when the incident photon energy is much less than the 
optical band gap, the Loudon-Fleury theory^^ has been 
widely employed to determine resonance profiles from 



2 



spin degrees of freedom via series expansions'"^, exact 
diagonalization of small clusters'^ or quantum Monte 
Carlo simulations'^^ of the Heisenberg model. Modifica- 
tions due to quantum fluctuations'^''^, bilayers'^, four- 
magnon processes''', couplings beyond nearest neighbor 
exchange'^, and ring exchange'^ have all been taken into 
account to give a thorough treatment of two magnon 
scattering from spin degrees of freedom in the nonres- 
onant regime. These approaches fail when the laser fre- 
quency is tuned to lie near an optical transition. In this 
regime, based on a spin-density-wave approach, Prenkel, 
Chubukov, and Morr have formulated a so called "triple- 
resonance" theory from which important features of the 
spectra can be derived'^. While good agreement was 
obtained for the resonant profile of the two magnon con- 
tribution to light scattering^, general features not related 
to the two-magnon peak are missed and lineshape calcu- 
lations are complex and only semiquantitative. 

An approach treating the full fermionic degrees of 
freedom is still lac;king. Recently exact diagonalization 
studies of the Hubbard model have been employed to 
yield lineshapes in the resonant limit from both spin and 
charge degrees of freedom". Yet the nonresonant and 
mixed terms were not taken into account. These calcu- 
lations also suffer the problems related to the finite size 
of the clusters (such as artificially broadening the delta 
functions to approximate thermodynamic-limit spectra). 
So generally, there is no theory for Raman scattering from 
both charge and spin degrees of freedom which predicts 
spectral lineshapes where both resonant, mixed, and non- 
resonant terms arc treated on an equal footing and do not 
suffer from finite-size effects. 

In this contribution we illustrate how to calculate the 
full electronic Raman response function, including con- 
tributions from the nonresonant, mixed, and resonant 
processes within a single-band model. Our model in- 
cludes interactions of the photon with all charge excita- 
tions of a correlated fermionic system, but does not take 
into account any scattering off of spin excitations. The 
scattering response is a complicated function of the cor- 
relations, the temperature, the incident photon energy, 
and the transfered energy. A short communication of 
this work has already appeared^°. 

Little is known about what the mixed Raman response 
looks like. We find that, as opposed to the nonresonant 
and resonant responses, which are manifestly positive, 
the mixed response is often negative (although the to- 
tal response always remains positive). The resonant re- 
sponse is expected to be large in the region where the 
transfered energy approaches the incident photon energy, 
called the double resonance, because the energy denom- 
inators of two pairs of the Green's functions in the bare 
response function approach zero. Interesting results are 
also anticipated in the strongly coupled (Mott-insulating) 
regime, when the incident photon energy is close to the 
interaction energy. Indeed, we find this is the case here. 
We also examine the situation where the initial photon 
energy is larger than the excitation energies in the cor- 



related band. This is the most common experimental 
situation in correlated materials with renormalized low- 
energy "bands" . The mixed and resonant responses also 
behave differently than the nonresonant response when 
we compare the Stokes (energy transfered from the pho- 
ton to the electrons) and the anti-Stokes (energy trans- 
fered from the electrons to the photons) responses. These 
are equal for nonresonant scattering, but the anti-Stokes 
response is much smaller than the Stokes response for 
the mixed and resonant scattering cases (introducing an 
asymmetry to the Raman scattering). 

The theoretical challenge in calculating the full inelas- 
tic light scattering response function is that the mixed di- 
agrams involve three-particle susceptibilities and the res- 
onant diagrams involve four-particle susceptibilities. It is 
only in the infinite-dimensional limit, where most of the 
many-particle vertex renormalizations vanish (all three- 
particle and four-particle vertices do not contribute; only 
the two-particle vertices enter) , can one imagine perform- 
ing the calculation of these susceptibilities exactly. It 
turns out that because the two-particle irreducible charge 
vertex is known exactly for the Falicov-Kimball modeP', 
one can calculate the full Raman response function in this 
case (since the general form of the charge vertex is not 
known for the Hubbard model, one can only perform ap- 
proximate calculations for that system even in infinite di- 
mensions; nevertheless, the diagrammatic analysis given 
in Section III holds for the Hubbard model, we just are 
not able to evaluate the expressions). 

We evaluate our exact expressions numerically and 
study their evolution as functions of the incident light 
energy and of the transfered energy. In the case of a 
correlated metal, we show how Fermi-liquid-like features 
evolve as the lifetime of putative quasiparticles decreases 
due to scattering. The results are even more interesting 
in the correlated insulator. We examine what happens 
to the isosbestic point identified in the nonresonant re- 
sponse, and how the presence of the charge gap affects 
the optical scattering. 

Inelastic light scattciring involves a coupling of photons 
to electronic charge excitations of the correlated material. 
The symmetry of the incident and scattered light relates 
to the symmetry of the charge excitations that are cou- 
pled to the light. There are typically three symmetries 
examined in experiments. The ^ig symmetry has the full 
symmetry of the lattice. This is measured, in a system 
with only nearest-neighbor hopping, by polarizing the in- 
cident and scattered light along the diagonal direction of 
the hypercubic lattice, so in large dimensions, we take the 
initial and final polarizations to be = = (1, 1, 1, ...). 
The Big symmetry is a d-wave-like symmetry that in- 
volves crossed polarizers along the diagonals. We take 

= (1,1,1,...) and = (-1,1,-1,1,...) for the Big 
channel. Finally, the i?2g symmetry is another rf-wave 
symmetry rotated by 45 degrees; it requires the polariza- 
tions to satisfy e' = (1,0, 1,0, ...) and = (0, 1,0, 1, ...). 
It turns out that the sector has contributions from 
nonresonant, mixed, and resonant Raman scattering, the 
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Big sector has contributions from nonrcsonant and reso- 
nant Raman scattering only, and the i?2g sector is purely 
resonant. This is generally true for a model on a bipar- 
tite lattice with nearest-neighbor hopping only. If longer- 
range hoppings are allowed, then all channels will have 
nonresonant, mixed and resonant contributions. 

While our approach towards analytic continuation is 
general, the overall complexity of the problem limits our 
evaluation of the light scattering cross section. By focus- 
ing on the spinless Falicov-Kimball model, we present a 
theory of light scattering from charge degrees of freedom 
only; valid for any incoming photon frequency. However, 
we are not able to address scattering from spin degrees 
of freedom resulting in two-magnon Raman scattering, 
for example. Nevertheless, we expect that our results 
help frame the physics related to resonance phenomena 
in paramagnetic correlated metals and insulators, and 
the behavior near a metal-insulator transition. 

Our plan of the paper is as follows. In Sec. II, we 
describe the general analytic-continuation formula that 
carries one from a time-ordered correlation function on 
the imaginary axis to the real response function. The 
formulas are completely general, and hold for the case 
of inelastic scattering of x-rays as well. The challenge is 
in evaluating the corresponding response functions along 
the real axis, which we know how to do only for the 
Falicov-Kimball model in infinite dimensions. In Sec. 
Ill, we evaluate the Raman scattering for the Falicov- 
Kimball model explicitly, calculating all response func- 
tions, and showing in detail how to perform all of the 
relevant renormalizations of the two-, three-, and four- 
particle correlation functions. In Sec. IV, we present our 
numerical results for Raman scattering at half filling. We 
examine the metallic case, the insulating case, and study 
the evolution of the Raman response as a function of the 
incident photon energy. We present our conclusions in 
Sec. V. 



II. GENERAL ANALYTIC CONTINUATION 
FORMALISM 

Our starting point is the expression for the inelas- 
tic light scattering cross section derived by Shastry and 
Shraiman^^ 



=27r^exp(-/3ei)<5(e/-ei-n) (2.1) 



for the scattering of electrons by photons of arbitrary 
wavelength (the repeated indices a and /3 are summed 
over) . Here = cj^ — ojj and q = ki kf are the trans- 
fered energy and momentum, respectively, while ijJi{f), 
ki(f^, and e^'^^^ denote the energy, momentum and polar- 
ization of the initial (final) states of the photons, £i(f) re- 
fer to the eigenstates describing the "electronic matter" , 



and g{q) = {hc^ /VLJqY/"^ is the "scattering strength" 
with Ljq = c\q\. Lastly, Z is the partition function. For 
an electronic system with nearest-neighbor hopping, the 
interaction with a weak external transverse electromag- 
netic field A is described by the following interacting 
Hamiltonian 



^-* = -^Ej'(^)-^(-^) (2-2) 

k 

g2 

+ 2fp:^2l^'^^-^^^°^'P^^ + k')Ap{-k'), 



where 



3a{q) = ^w«(fe)4(fe + ql2)ca{k - q/2), 



de{k) 
dka 



(2.3) 



are the current operator and Fermi velocity, respectively, 
and 



= E + 9/2)^-(fc - 9/2) (2-4) 



is the so-called stress tensor. As a result, the scatter- 
ing operator M{q) has both nonresonant and resonant 
contributions 



f M'^'^iq) i)={f\'y^M\i) (2.5) 



+ 



£l- £i- LOi 

{f\u-km {I \30{kf)\z) 

El -Ei+OJf 



with the sum I over intermediate states. The term with 
the stress tensor is the nonresonant contribution, while 
the term with the square of the current operator is the 
resonant contribution. Now the Raman-scattering cross 
section contains nonresonant, mixed, and resonant con- 
tributions (because it is constructed from the square of 
the scattering operator): 

R{q, n) = RN{q, n) + RMiq, O) + RRiq, O), (2.6) 

where the nonresonant contribution is 



RN{q,n)^2TTg^ki)g^{kf) 
exp(-/3£j) _ 



(2.7) 



%f 7/,i H^f -£i-^), 
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the mixed contribution is 
RM{q,^) = 27rg^iki)g^ikf)Y^ 

i,f,l 



jf,i Ji,i 



cxp(-/3£i) 



Jf,lJl,i 



(2.8) 



ei - Ei - u)i + iO+ si - £i +LOf - iO+ 



+ 



Ji,l JlJ 



+ 



Af)Ai) 

Ji,l JlJ 



ei - Ei - Ui - iO+ si- 6i+u>f + iO+ 

x(5(£/ -Si- n), 



lf,i 



and the resonant contribution is 

Ruiq,n) = 2ng\h)9\kf) ^-^^^^ (2.9) 

i,f,l,l' 

( Ai)Af) Af)A^) \ 

^ [ JlJ |_ Ji,l ^IJ \ 



ei — Ei — LOi — iO+ El — Ei + Uf + ?0+ 

Af)Ai) A^ Af) 

Jj\i'Ji',i Jf,i'Ji',i 



\£l'-£i-LOi+iO+ El' — Ei + U!f - iO+ I 
X 6{Ef - Ei - fl). 

In these equations, we have introduced the following sym- 
bols 



7 = Xle^7«,/3(q')e^, 



a/3 



(2.10) 



^(^iMkf), 



\A\ f) for the matrix elements 



with the notation Aij 
of an operator A. 

In general, the matrix elements that enter into 
Eqs. (2.7-2.9) are not easy to calculate for an interact- 
ing system, so the summations are problematic to eval- 
uate. Instead, these expressions usually are evaluated 
via Green's function techniques starting from correlation 
functions evaluated on the imaginary axis and then per- 
forming an analytic continuation to the real axis to get 
the physical response functions. This procedure becomes 
more complicated when the number of matrix elements 
that enter into each term in the summations increases, 
because it requires the evaluation of a more complicated 
correlation function on the imaginary axis. Our strategy 
is to first consider the analytic continuation procedure in 
a general sense, which holds for any model Hamiltonian 
and for arbitrary momentum transfer. We will derive 
connection formulas between the Matsubara frequency 
axis correlation functions and the analytically continued 
response functions on the real axis. But those expres- 
sions will require us to be able to evaluate a number of 
different susceptibilities, and those expressions are not 
known for arbitrary Hamiltonians. We will show how to 



evaluate them exactly for the Falicov-Kimball model in 
the next section. 

We derive below the connection formulas between the 
imaginary-time response functions and the real frequency 
response functions for the general case. In all of our nu- 
merical results, we examine only optical light scattering 
where we approximate fej =kf = q = 0. 



A. Nonresonant scattering 

The nonresonant scattering in Eq. (2.7) is proportional 
to the spectral density function. The spectral density 
cannot be calculated directly but is instead obtained from 
the analytic continuation of the imaginary-time response 
function constructed from the time-ordered product of 
two stress-tensor operators 



xi'Ur,r') = (r.7(r)7(r')> 



^7.7 



(2.11) 



with the T dependence of the operator determined by the 
Hamiltonian in the absence of the electromagnetic field 
(the symbol 7^ denotes time ordering). The first step 
is to calculate the double Fourier transformation to the 
Matsubara frequency axis 



/3 /3 




(2.12) 

for Bosonic Matsubara frequencies = mT2n with 
(3 — 1/T. In thermal equilibrium, the two-particle corre- 
lation function depends only on the difference of the two 
time variables, so the double Fourier transform becomes 
a "diagonal" function, evaluated as 



ij 



exp(-/3£^) 7»,/7/,i 

Z £f — Ei — if 



(2.13) 



X [1 - exp(/3(£i - £/))] . 



In order to extract the spectral density of states from the 
Matsubara correlation function in Eq. (2.13), we perform 
the analytic continuation ii/ — > O ± iO+ which yields for 
the nonresonant scattering the known expression 

where we introduced the nonresonant response function 

r{q, = ^ {x'-'M-^ - *0+, n + *0+) (2.15) 



-x^!^(-fi + io+,n-io+)} 



evaluated on the real axis. A similar strategy is used 
to determine the mixed and resonant contributions as 
described in the next two subsections. 
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B. Mixed scattering 

In the case of mixed scattering in Eq. (2.8), the cal- 
culation begins with the multi-time correlation function 
constructed from the stress tensor and two current oper- 
ators 

= (^r7(r)j(^HT')j^'Hr")) • (2.16) 

We define the Fourier transform as before, with respect 
to three Matsubara frequencies (all with the same sign 
of the exponent). Once again, in thermal equilibrium we 
have imaginary-time-translation invariance, so the sum 
of the three Matsubara frequencies must vanish, yielding 



1 



Z 

exp(— /3£i) 



(2.17) 



{sf -£i + wi){ei - £i- ius) 

exp(-/3£/) 



{ei -ef + iv2){ei - £/ - iyi) 

exp(-/3e,) 



hf,i 



{si -ei+ iv3.){ef - £i - W2) 
exp(— /3£j) 
(£/ - Si- ivi){ei - ei + ivs) 

exp(-/3£/) 



{ei - ef - iv2){£i - £/ + ivi) 
exp(-/3£() 



(£i -£i- iv2,){sf - ei+ W2) 



which contains 3! = 6 terms collected into two groups of 
terms connected by cyclic permutations, with 



(2.18) 

After analytic continuation iu^ Za with the constraint 



.(3) 



Zl+ Z2 + Zz = 0, 



(2.19) 



one can sec that the expression in Eq. (2.17) has three 
branch cuts when Imz^ 0^ (for a = 1, 2, or 3). Note 
that the constraint in (2.19) forbids only two of the z^s 
to simultaneously have Im^^ = 0^, but the imaginary 
part of all three can vanish simultaneously. In order to 
produce the expression for the mixed Raman cross sec- 
tion in Eq. (2.8), we need to focus on the branch cuts that 
occur when z\ — ± zO+ and Zi ^ Sl± iO+ in order to 
produce the appropriate (S-function and matrix elements 
in the mixed scattering cross section. The correspond- 
ing discontinuity across the branch cut when Imzi = 0, 
occurs when the terms in Eq. (2.17) are analytically con- 
tinued with zi moving onto the real axis. In the first case 



when zi — ft ± iO+, and 2:2 — > — Z3, we find 

zi— >-n-io+ 



1 



2m 



(2.20) 



1 

27ri 
/3) 



X^'},,(-r!-z0+,r!-Z3,^3) 



= (l-e-^")E 

X (5(£/ - £i - O) 



Jf,l Jl,i 



+ 



Ai)Af) 
Jf,iJi,i 



Sl-Si-Z3 ei-Ef + Zs 



and in the second case when zi ^ fl± iO+, and Z2 
—fl — zs, we find 



-xfj,i{zi,Z2,Z3) 



27ri 



zi^n-io+ 



(2.21) 



xfj/^ + lO^,-^-Z3,Z3) 



2-Ki 

.(3) 



- xi,^,,(fJ-iO+,-f2-03,^3) 



Jf,l Jl,i 



+ 



Ai)Af) 

Jf,lJl,i 



Si - Si- Z3 ei-£f + Z3 



X 5{£f -£i + O). 



The sum on the right-hand side of Eq. (2.20), with 
Z3 = Ui — iO"*", is proportional to the first two terms in 
Eq. (2.8). The sum on the right-hand side of Eq. (2.21), 
with an interchange of z / in the summation and 
Z3 = —LUi — zO+, is proportional to the last two terms 
in (2.8). Hence, we arrive at the general expression for 
the mixed scattering 



RM{q,^i) = / XM(g,") 



(2.22) 



1 - exp(-/31^) 
with the mixed Raman response function defined by 

(2.23) 

- X~%i-^ + *0+, -ujf + iO+,u;^ - iO+) 
+ X~ j,i(f^ + iO+,LUf + 1O+, -w, - zO+) 

- xff,i{ft -iO+,iJf+iO+, -Ui - iO+) 

on the real axis. The operators 7(9), j^.''^(fe/) andj^*^(fei) 
are Hermitian for optical light scattering, which has van- 
ishing momentum ki=kf = q = 0, and, with the use of 
Eq. (2.18), the expression in Eq. (2.23) can be rewritten 
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as: 



with 



XM(fi) =XM{q = o,n) 



(2.24) 



2m 



xfj,i{^ + LOf - iO+, -Ui + iO+) 



+ xfj,i{^ + ^0+, w/ + iO+, -iJi - iO+) 

(3) 



C. Resonant scattering 

For the resonant scattering case in Eq. (2.9), the pro- 
cedure is similar: one has to calculate the multi-time 
correlation function constructed from the four current 

operators 

X^lu{ri,r2, T3, T4)= {Trj('^ (ri)j(^) (r2)j(/) (r3)j(^) (r,) 

(2.25) 

Once again, defining the Fourier transform in terms of 
four Matsubara frequencies (with the same sign in the 
exponent) yields the following result (with the delta func- 
tion arising from the time-translation invariance) 

Xijj,Si^i, i^2, ii^3, i^i) =5{ux+U2 + u^ + u^) (2.26) 

xljj,ii'ii^i,i'^2, ifi, ifi) + x[j^ij-{ii^i,ii^2,ii^i, i^s) 
+ xtlfjii^i^i^i^iv^^ivs) ■ 
Hero wc introduce the generic four-particle susceptibility 

XA,B,C,£>(*^l'^^2,iz^3,«Z^4) = ^ Ai^iBijCf^i'Di^^i^ 

i,f,l,l' 

(2.27) 

cxp(-/3£i) 

{ei -ei + wi){ei' - Si- W4){£f - - - wa) 

exp(-/3£;) 

(e/ - si + iv2){£i - £i- - £i- ivi - ifi) 

exp(-/3£/) 

{ei> - £/ -h W3)(e; - £/ - W2){,Si - £/ - ii'i - i^2) 
exp(-/3£;/) 



-I- 



+ 



+ 



[si - Si' + ivi){ef - £i' - W3){ef - sv - 1^2 - _ 



exp(-/3£i) 



-I- 



-I- 



+ 



(£; - £i - iVl){£l' - £i + iV4.){£f - £i + + Wa) 

exp(-/3£;) 

(£/ - £; - iu2){£i - ei+ iy\){£i' - £1 + iv^ + ivi) 

exp(-/3£/) 

{si' - £f - if3)(£; - £f + W2){£i - £/ + Wi +^2) 
exp(-/3£;/) 



{£i - £v - iVi){£f - £l' + ii^s){£f - £i' + ii^2 + ii^s) 



Xa,b,c,d(«^1' «^2, it's, ivi) 



(2.28) 



The expression in Eq. (2.26) contains 4! = 24 terms col- 
lected into six different groups of the terms, with each 
group member connected by the cyclic permutation of 
four objects. 

After analytic continuation iua — > Za with the con- 
straint 



Zi + Z2 + Z3 + Z4 = 0, 



(2.29) 



one can see that the expression in Eq. (2.26) has branch 
cuts when any luiZa 0^ or when any pair lm.{za + 
zp) — > 0^. The (5-function in the expression for the reso- 
nant scattering cross section in Eq. (2.9) is connected to 
the branch cut at Z3 -|- 04 = — zi — Z2 — > ± iO+ and the 
discontinuity of the response function across this branch 
cut is equal to 



1 



l^.xf}j,i{zuZ2,Zz,ZA] 



2m 



z^-\- Z4 — — zi— Z2—*^~i—iO+ 



(2.30) 



(1 



-/30 



)E 



-6i£f- 

Ai)Af)Af)Ai) 

JlJ Jf^l'-IVi 



{£l - £i + Zl){£l' - £i - fl + Z3) 



+ 



Af)Ai)Ai) AS) 

Ji,l Jl,fJf,l'Jl',i 



(£;-£,- n - Zi)(£r - £i - Z3) 



Ai)Af)A^) Af) 

Ji,l JlJ Jf,l'Jl',i 



+ 



{£1 -£i + Zi){£i, 

AS)A^AJ)Ai) 

Ji,l JLfJf,l'Jl',i 



Z3) 



{£l-£i-n- Zi){£i> 



Z3) 



The analytic continuation procedure then requires us to 
take the following limits 



zi — 

Z2 - 
Z3 - 
ZA — 

and then take the limit 



■ -Wj - iO+, 

■ LUf + iO+, 
-oj'f + iO+, 
J, - iO+ 



(2.31) 







(2.32) 



in order to reproduce an expression proportional to the 
resonant scattering cross section in Eq. (2.9). The final 
general expression for the resonant scattering becomes 

with the resonant Raman response function defined by 
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(4) / \ 



Z3 + Z4 = — Zl—Z2^O-\-i0'^ ' 



Z3+Z4 = — Zl— 22— ^f^— iO+ 



zi- 

Z2- 
Z3- 
Zi- 



■-Wi-iO+ 

-U!'j,+i0+ 
^aj'-iO+ 



(2.34) 



jtJ^—LOi — ♦O 



note that it is critical to perform the analytic continua- 
tion of Z3 + Z4 = —Zi — Z2 ^ iO+ first and then ana- 
lytically continue the other frequencies [as in Eq. (2.30) 
and (2.31)] since these procedures do not commute with 
one another. 



III. EXACT RESULTS FOR THE 
FALICOV-KIMBALL MODEL 

We now evaluate the general expressions derived above 
for the case of optical Raman scattering, where all mo- 
menta vanish (fej = kf = q = 0) and for the spinless 
Falicov-Kimball model. The Falicov-Kimball model in- 
volves the interaction of conduction electrons with local- 
ized electrons and has the following Hamiltonian^-'^ 



with 



{i,3) 



~ 1^ ^i4^i +Wi) + U^ c\ci 



Wi 



(3.1) 



where cj (cj) create (destroy) a conduction electron at 
site i, tVi is a classical variable (representing the localized 
electron number at site i) that equals or 1, is a renor- 
malized hopping matrix that is nonzero between nearest 
neighbors on a hypercubic lattice in D-dimensions (and 
we take the limit f — » oo), and U is the local screened 
Coulomb interaction between conduction and localized 
electrons. denotes a sum over sites i and nearest 

neighbors j. Ef and fi are adjusted to set the average 
filling of conduction and localized electrons. In our cal- 
culations the average filling for each is set to 1/2, respec- 
tively in = U/2, Ef = 0). 

This model can be solved exactly by using DMFT, as 
first described by Brandt and Mielsch^'^. The algorithm 
used to solve for the local Green's function at site i, de- 
fined by 



Gi{T) = -{TMT)clm 



(3.2) 



where the angle brackets denote the trace weighted by 
the Boltzmann factor exp[—/3H]/Z. We usually work 
with the Fourier transform of the imaginary-time Green's 
function to yield the Matsubara frequency Green's func- 
tion. The momentum-dependent Green's function be- 
comes 



1 



efe 



(3.3) 



D 



efe 



= — lim — ;= cos ka 



(3.4) 



* a=l 

being the noninteracting band energy, and 

Zm = i^m +1^- ^m- (3.5) 

The local self-energy is a solution of the following set 
of equations: 



Gr, 



1 ^ 1 _ 1 



(3.6) 



Wi 



+ ■ 



1 — Wi 



iu>m + H — \m — U iWrn + M ~ -^m 

where we introduced the self-consistent dynamical mean- 
field of Brandt and Mielsch (denoted A); the self-energy 
can be expressed as a simple function of this field 



S„ = Uwi + - 



iu>m -\- jJL- \m-U{l- Wi) ' 



(3.7) 



Here wi is given by Wi = ^^^f '^^^2o{U — ij)/Z, with 



A. Nonresonant scattering 




FIG. 1: Feynman diagrams for nonresonant Raman scatter- 
ing. The wavy lines denote photon propagators and the solid 
lines denote electron propagators. The cross-hatched rectan- 
gle is the reducible charge vertex. In the -Big channel, only the 
bare (first) diagram enters, while in the Aig channel both di- 
agrams enter. The symbol 7 denotes the stress-tensor vertex 
of the corresponding electron-photon interaction. 

The case of nonresonant Raman scattering was con- 
sidered by two of US'", so we only sketch the derivation 
to show our notation and to present the final results. In 
general, the two-time correlation function in Eq. (2.13), 
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constructed from the stress operators, can be represented 
by the summation of a generalized "polarization" 

X^Jyi-i^yi^) = T^U{i(Jm.,i(^m+,^) , (3.8) 



where we use a shorthand notation iu, 



iuJrn + iv 



and Gm+v = G{iLOm+v) = G{iLOm + w) and similarly for 
S and Z. A tedious calculation shows that'' 



in the B-^g channel and 



2 Zi^j^u -^m 



(3.9) 



1 ^m+L' 



(3.10) 

in the A^g channel (see Fig. 1 for the relevant Feynman 
diagrams) . 

Since one can show that there are no additional singu- 
larities or non-analyticities in Eqs. (3.9) and (3.10) con- 
nected with the denominators,^** one can directly perform 
the analytic continuation and replace the sum over Mat- 
subara frequencies in Eq. (3.8) by an integral over the 
real axis 



+ 00 



(3.11) 



X [n(a; - iO+, w + iiy) - n(w + iO+, w + ii^) 
+ n(w -iu,ui- iO~^) - II{u -iv,LO + iO+)] , 

where /(w) = 1/[1 + exp(/3a;)] is the Fermi distribution 
function. After substituting Eq. (3.11) into the expres- 
sion for the non-resonant response function in Eq. (2.15), 
we obtain 



+00 

Xiv(l^) = ^^ / Mf{^)- f{^ + m (3.12) 

— oo 

X Re {n(w - iO+, UJ + Q. + zO+) 
-n(a;-iO+,a;-|-fi-iO+)}. 



Now we can take the trivial analytic continuation of 
Eqs. (3.9) and (3.10) to find the final expressions for the 
nonresonant Raman response function: 



in the B\g channel, and 

+ 00 

XN,A.,{^) = ^ j du[f{^)-f{^ + m (3.14) 



X Re 



[■Ejco) -T,*{u) + Q)][Z*{u} + O) - Zjco)] 

G-'^*{uj + n)-G-^uj) 
[£{oj) - E{oj + n)][Z{uj + n)- Z{oj)] 



G-i(w + 0)-G-i(w) 
in the Aig channel, respectively. 

B. Mixed scattering 



The mixed Raman response corresponds to the scat- 
tering processes that involve three external vertices: one 
stress tensor and two current operators, and there are 
two types of diagrams corresponding to the direct and ex- 
change processes (see Fig. 2). There is no mixed Raman 
response for the channel because the stress tensor 
vanishes for the case of nearest neighbor hopping only. 
In the Big channel it appears to be only a bare response 
(we will sec below that it actually vanishes) and for the 
Aig channel, the bare mixed response is renormalized by 
the irreducible charge vertex. 



+ T 




FIG. 2: Feynman diagrams for the mixed contributions to Ra- 
man scattering. The symbols j/ and ji remind us to include 
the relevant vertex factors from the current operator in the 
electron-photon interaction. The mixed contribution vanishes 
in the B2g channel, it consists of only the bare diagrams on 
the top line in the Big channel (and turns out to be a 1/D 
correction), and all diagrams enter for the Aig channel. 



2 

XN,B,,{fl) = ^ J dcv[f{uj)-f{u; + n)] (3.13) 



1. Big channel 
il, the mixed R; 

Z*{ui + — Z{u!) Z {oj + — Z (uj) j tains only the bare direct and exchange contributions 



^ , Giu;)-G*{u; + n) _ G{uj) - G{uj + O) , channel, the mixed Raman response con 
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(first two terms in Fig. 2, respectively): 



2. Aig channel 



m k 

D D D 

X ^(-1)" COS fee y^(-l)^ sinfc/3 sin fc-y 

a=l 13=1 7=1 

(fc) 

(fe)], 



because the symmetry of all two and three-particle ver- 
tices is that of the lattice (^ig), so all rcnormaliza- 
tions vanish (recall the current operator has odd parity, 
whereas Cfe is even in k). The expression in Eq. (3.15) 
has nonzero values only when the subscripts are equal 
a = /? = 7. In this case, wc expand the product of 
Green's functions into partial fractions over and the 
summations over momentum involve only expressions of 
the type [Cm = -sgn(ImZm)]: 



I t*^ " 1 

h a=l 



(3.16) 







'It 1 i cos h 

— sm k cos k e ^ 
27r 



with Jo being Bessel's function. The last exponent is 
expanded in a power series over that yields in the 
D ^ OQ limit 



*4 



lim — — — 



C™^e-^^/**%rfc(.C 



t* 



t* 



j.*2 



= n^^"^ 175 - 1) - 

D— >(x) 4i/ 



SO the mixed contribution vanishes in the Big channel. 



In the Aig channel, the mixed Raman response con- 
tains both bare and renormalized contributions: 



(3.18) 



f 1 t*2 ^ ^ 

m \ k a=l /3=1 

(fe) 

1 ^2 ^5 

(fe) 

xTf {iujm, ii^m + ivi - iff) 

^■L^^kGm{k)Gm+Ui-ui{k) \ 



N 

+ 



k 

iUf -iVf 



The renormalizations are only with respect to two- 
particle vertices, because the current operators are odd 
in parity and cannot be renormalized by a local three- 
particle vertex (note that wc cannot provide a general 
proof that the relevant three-particle vertex is local, but 
a strong-coupling analysis indicates this is so to lowest 
order). Here 



f {iuim, iijJm + m - iff) 



(3.19) 



1 — rr (iuJrn, i^JJ' 



is the total (reducible) charge vertex. In the D — 00 
Falicov-Kimball model, the irreducible charge vertex sat- 
isfies 



TGr, 



(3.20) 



Gm+iy 

on the Matsubara frequency axis^'^. Substituting into the 
expression for the reducible charge vertex gives 



r {iu}m,i^m+Vi-Vf) 



1 Zm+vi 



-'m+L'i—Uf 



T [ G 



Zm, — Z. 



m+Vi — i/f 



m+Vi — Vf 

(3.21) 
2 \ 



Gm+Vi-Vf {il^i — il'f){Gm — Gm+Vi-Vf) 



Now Eq. (3.18) has nonzero values only when a = (3 

and, noting that in the D ^ 00 limit one can replace 
(3.17) 

sin ka by its average value ^, yields 

x\j^iiin-ivf,ii'f,-ivi) (3.22) 

=T^ |^n*^^^(«Wm - Wf,iLOm + iVi - ivf,iuim) 
m 
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where 



Here the analj^tic continuation of Eq. (3.23) is 



(3.23) 



*2 



-'m+Vi — Vf 



2{iVi — iff) Gm — ^m+Ui-Uf 



G 



m+Vi — Vf " 



Grt 

'm—i/f 



G„ 



G 



ra+Ui — i/f 



Gri 



-'m—Vf 



-'rri—Vf 



-'Ta-\-Vi—Vj 



In the case when there are neither singularities nor non- 
analyticities in Eq. (3.23) connected with the denomi- 
nators, one can trivially perform the analytic continua- 
tion and replace the sum over Matsubara frequencies in 
Eq. (3.22) by an integral over the real axis: 



+00 



n^'^^ (w — ivf, bj -\-ivi — ivf, uj — iO"'") 

{oj — iuf, ui + iui — ivf, ui + iO'^) 
{w — iui,u} — iO'^,u} — ivi + iuf) 
{ui — iui,u} + iO'^,u} — ivi + iff) 
{lj + iuf,ui — iui + ivfjUi — iO'^) 
{u! + iuf,ui — ivi + ivfjUi + iO"*") 
(w + if,, w — iO'^,u} + ivi — iff) 
{co + iui^LO + iO"'", w + ivi — iff) 
(w — iO'^,u) + ivi, oj + ivf) 
(w + iO"*", w + if j, LO + ivf) 



_n(3 

_n(3 
- n(3 



(cij — iO"'', uj ■ 



iUf) 



{co + iO'^,LO — ivi, oj — ivf) 



Only the first eight terms contain the difference ivi — ivf 

and hence contribute to the mixed scattering. Substitut- 
ing Eq. (3.24) into Eq. (2.24) we get the final expression: 



-t-00 



25) 



XM,A.,(1^) = / Mf{^)-f{^ + m (3- 

— CSO 

X Re |n(3) (uj - Ljf + iO+ ,uj + n + iO+ ,10 - iO+) 

- n(3) {uj - ojf + iO+ ,uj + n - iO+ ,uj - iO+) 
+ n(3^ (w - w/ - iO+, w + O + iO+, UJ - iO+) 

- n^^^ {uj-uf- iO+, w + O - iO+, UJ - iO+) 
+ n(^) {uj + uJi-iO+,uj + n + iO+, UJ - iO+) 

- n(^) {uj + uJi-iO+,uj + n- iO+, UJ - iO+) 
+ n(^) {uj + uJi + iO+,uj + n + iO+, UJ - iO+) 

- n(^)(a; + a;i + iO+,a; + n-iO+,a;-iO+)|. 



n^^) {UJ1,UJ2,UJ3) 
+*2 \ 



S(u^3) 

Us) G{U3) 



2{UJ2 

- G{^^js)-G{}^ 



- S(C^2) 
G{UJ2) 



(3.26) 



G(a>3) - G(c^i) . G{uj2) - Gjuj,) 

-ZiuJ3) ^ Z{UJI)-Z{UJ2), 



C. Resonant scattering 

The resonant Raman response corresponds to scat- 
tering processes that involve four external current ver- 
tices. The correlation function constructed from four 
current operators contains six types of diagrams corre- 
sponding to the different direct and exchange processes 
(see Fig. 3). It should be noted that since j*^*' and j'^^^ 
are odd functions of momentum, the only way to get a 
nonzero momentum summation is to have an even num- 
ber of current operators in any given momentum integra- 
tion ("current-operator pairing"). Hence all local three- 
particle and four-particle vertex renormalizations must 
vanish, although two-particle vertex renormalizations are 
possible. For the Big and B2g channels the "current- 
operator pairing" is possible only between either both 
incoming j^^^ or both final j^^^ current operators, but 
for the Aig channel all operators can be involved in the 
"pairing" and the contribution from the bare diagrams 
in the Aig channel is three times larger than for the Big 
channel. As a result, in the Big and B2g channels we have 
contributions from the first two diagrams in the first four 
lines and from only the first diagram in the last two lines 
of Fig. 3, and in the Aig channel all diagrams contribute. 

For the Big and i?2g channels, the product 
is orthogonal to the charge vertex with Aig symmetry, 
so the diagrams are not renormalized across the vertices 
that contain both j^^^ and j^^^ factors. In addition, for 
the B2g channel, the polarization vectors select either odd 
or even momentum coordinates and, as a result, the res- 
onant Raman response for the i?2g channel is four times 
smaller than for the Big one, and it is the only contri- 
bution to the total Raman response in the B2g chan- 
nel. In the Aig channel, besides the diagrams presented 
in Fig. 3 that include all possible horizontal and verti- 
cal "ladder" renormalizations, one could renormalize by 
parquet-like terms that involve simultaneous horizontal 
and vertical renormalizations. But it can be shown (see 
the Appendix), that such contributions are 1/D correc- 
tions, and disappear in the — > 00 limit. 

As a result, the Fourier transform of the four-time cor- 
relation function constructed from the current operators 
can be represented in the following form (ifj — ivf = 



11 



AA, 



ilililll 



.X, .X 



+ 



X" Y X 

>f X 



- T 




- T 



Jf Jf 



\ 



Jf Jf 




FIG. 3: Feynman diagrams for the resonant contributions to 
Raman scattering. Only the first two diagrams in the first four 
lines, and the first diagram in the last two lines contribute 
in the Big and sectors. The response includes all 
diagrams. 



X}jj,i{-m,i'^f,-if'f,ii'l) = 



(3.27) 



mi 



In the Big and B^g channels, n^^-* contains only the 
bare contribution (corresponding to only the first dia- 
gram on last two lines of Fig. 3) 



(3.28) 



(4) 



where 



nbtre(^'^l'*'^2,«W3,ia;4) 



(3.29) 



D D 
a=l p=l 



'D^N^ (Zi - Ck){Z2 - ek){Z3 - efe)(Z4 - Ck) 



j.*4 " 



iZ2 - Zi){Z3 - Zi){Zi - Zi) 

G2 



{Z^- Z2){Z3- Z2){Z^- Z2) 

G3 

(Zi - Z3){Z2 - Z3){Z4 - Z3) 



[Zi - Zi){Z2 - Zi){Z3 - Zi) 



However, the other polarization 11^^' contains a vertical 
"ladder" renormalization (corresponding to the first two 
diagrams on the first four lines of Fig. 3) 



\].^pg^^{iijji,iw2,iu:3,iu:i) (3.30) 
= 'n.^l.^{iuJi,iuJ2,ii^3,i<^A) + 'n.[^\iuJi,iu)2,iu)3,iu>A), 
n/%2g(«wi,ia;2,«W3,iw4) = ^^f^B^^{iL0i,iuj2,i'^3,ii^i) 
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with 



(3.31) 



= - ^ ( ^ E ^ E Si^' fcaGl(fc)G2(fc)G4(fc) J 
V fe a=l / 

xf (2, 4) ( ^ ^ ^ ^ sin2 fc„G2(fe)G4(fc)G3(fc) 



a=l 



^t^^G4-Gi G2-Gl^^ r(2,4) 



A Zi - Z2 ; (^2 - ^4)2 

G4 — G3 G2 — G3 

■^3 — ^^4 -^3 ~ ^^2 



Using the solution of Eq. (3.21) in the Bethe-Salpeter-like 
equation (3.19) yields 

'n!^^ {iL0\,iw2,iuJz,iu)i) = n^^^(iwi,iw2,iw3,iw4) (3.32) 

^fr^^ [iuji , iw2, «W3, ia;4) 



+ 



IUJ2 — lt04 



with 



n*;''-' {ioji, iu}2, ioJa, ioja) 



x*4 



4 (G2-G4)(^2-^4) 



(3.33) 



G/^ \ / 



Z\ — Zi — Z2 J \ Zs — Z4, Zs — Z2 J ^ 



and 



t 



*4 



1 



4 G2-G4 

G/^ \ / 



(3.34) 



Z\ — Z/^ Zi — Z2 J \ Z^ — Zn Zy, — Z2 

^'^^ (iijJ^^iL02-,iiO\^iu}A) = — \l/*;''^(iwi, ZW4, »W3, iW2). 



In the A\g channel we have contributions from all the 
diagrams in Fig. 3, hence 



(za;i,ia;2,«a;3,za;4) = n^j^^^^ {iuJi,i(J2,i(^3,ioj4,) 



(3.35) 

= 3n[^^g (i wi , ia;2 , «W3 , i W4 ) + n^^^iwi , ia;2 , « W3 , ia;4 ) 

Here the last term corresponds to the horizontal "ladder" 
renormalization (the last diagram on each line of Fig. 3). 

Next, we perform the analytic continuation in 
Eq. (3.27) and replace the sum over Matsubara frequen- 
cies by an integral over the real axis in the same way as 
was done in Eq. (3.24) for the mixed scattering. Then 
we substitute it into the expression in Eq. (2.34) for the 
resonant Raman response. After some tedious algebra, 
we achieve the final expression for the resonant Raman 
response of the D = 00 Falicov-Kimball model: 



+00 

^^^^^ = (2^ / '^^^^''^ ~ ■^^'^ + ^^-^^^ 

— 00 



{nf(-- 


- iO+,u- 


-UJf 


-zO+, 


LU - 




\- iO^ , LO — LUf - 


f iO+) 




- iO+,w- 


-U)f 


-iO+ 




+ n 


— iO^ ^u) — LOf 


+ iO+) 


+ nf\uj 


- iO+,uj- 


f iOi 


-iO+, 


U! 


+ n + iO+,Lu + iUi- 


f iO+) 


-nr'(. 


- iO+,w- 


f LOi 


-iO+, 


LU 


+ n 


— iG^ ,LU -\- CUi - 


f ?0+) 


+ ufX- 


- iO+,uj- 


f 


+ iO+, 


LU 


+ n- 


+ iO'^,Lu — UJf 


+ iO+) 




- iO+,uj- 


f LOi 


+ ^0+, 


LU 


+ n- 


— iO~^,LU — UJf 


+ ^0+) 


+ ufX^ 


-iO+,w- 




-iO+, 


LU 


+ n- 


+ iO~^,LJ — Wf 


-zO+) 


- li^h- 


— iO^,LO + Ui 


-iO+ 




+ n 


— iO~^,cu — u)f 


-iO+) 



} 



The analytic continuation in Eq. (3.36) can be found simply by substituting itUa ^ ± «0+ in the corresponding 

expressions in Eqs. (3.28)-(3.35) which will not be explicitly repeated here. It might appear that the first four terms in 
braces contain divergences connected with vanishing denominators in the last term in Eq. (3.32), but the contribution 
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of these terms into the expression in braces in Eq. (3.36) must be considered in the limit: 

hm — — iO+,a; + n + iO+,a;-a;/- A + iO+) (3.37) 

- '^['^^ (w + iO+, w - w/ - iO+, w + 17 + iO+, w - w/ - A + iO+) 

- ^-^^^ (w - iO+, w - w/ - iO+, w + O - iO+, w - w/ - A + iO+) 
+ (w + iO+, w - w/ - iO+ , w + O - iO+, w - w/ - A + iO+) 
+ (w - iO+, w + + A - iO+, w + f2 + zO+, w + Wj + iO+) 

- (w + iO+, w + + A - iO+, w + + iO+, w + + iO+) 

- (w - iO+, w + Wj + A - iO+, w + n - iO+, a; + + iO+) 

+ *^^^(a; + iO+, w + + A - iO+, w + - zO+, w + + iO+) | , 

where A = uj'^ — ujf = ui[ — uji. When the hmit A ^ is taken, we find that the imaginary part of Eq. (3.37) 
diverges, but the real part (which is all that contributes to the Raman scattering) is finite and can be calculated using 
I'Hopital's rule: 



-iRe{^'C 



with 



and 



uj - iO+, Lo~LOf- iO+, w + O + iO+, Lo-ujf + iO+) (3.38) 

- ^'{uj - iO+,a; - w/ + iO+,w + - iO+,a; - w/ - iO+) 

- *'(a; - iO+, w - - iO+, a; + 0-iO+,a;-w/ + iO+) 
+ *'(a; - iO+, a; - w/ + iO+, a; + + iO+,a;-a;/ - iO+) 
+ *'(a; - iO+, w + + zO+, w + O + iO+, w + - iO+) 

- *'(a; - iO+, w + - zO+, w + O - zO+, w + + iO+) 

- *'(a; - iO+, a; + + iO+, w + O - iO+, w + - iO+) 

+ ^'(w-iO+jW + Wj -iO+,a; + + iO+,w + Wi + iO+)| , 

I 



f(M,3,4) = f?iqMM|l (3.39) f = l-f. (3.42) 

When the ground state of the "paramagnetic" phase is 
metallic, the derivative of the self energy is straightfor- 
d^f'{l, 2, 3, 4) _ dGi/dZi ward to calculate as 

dZ^ ~ ~T (G2 - G^Y rfS_ 2 1-ZG 

(3.40) i*2 G2 ^ ' 
Gi-Gi G2-Gi \ (G^-Ch G2-(h\ X {T.-Uwrf 

_ 1 / dGj/dZj ^ G4-G1 \ but in the "insulating" phase, when U'^Wi (1-Wi) > i*^/2 

4 G2 — G4, \ Zi — Z/^ {Zi — ZiY ) filling one has to include the contribution from 

' G — G G ~ G \ delta-function peak in the imaginary part of the self- 

energy which yields the additional contribution 



Z2, — Zi Z-i — Z; 



t*^ 1 f G^-G^ G2-GA ^ = 2 1|G 

' ' i*2 G2 ^ ^ 



4 G2 — G4 y Zi — Z4 Z\ — Z\ 



(S-t/wi)2 

( dG^/dZi , g4-G3 ^ X _ + (s _ C7«;i)'(l + F^^) 



Z2, — Zi {Z3 — Zi) 



dG 2 

= ~ ■^^)' (3-41) to the derivative. 



±m [U'^wi{l-wi)-t*'^/2] 



dSjuj - U[l - wi]) 
duj ' 
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D. Bare contributions and multiple resonances 

In summary, the total Raman response function is the 
sum of the nonrcsonant [Eq. (3.12)], mixed [Eq. (3.25)], 
and resonant [Eq. (3.36)] contributions and has a com- 
plicated form. It is educational to consider the contribu- 
tions of the bare diagrams, which can be summed up and 
rewritten in the following form^^: 

^(^) = ^E / div[f{u,)-f{u, + n)]Ak{iv)Ak{iv + n) 



(3.45) 



where 7fc 



^fc(w) = ^ ImGfc(a; — zO+), and 
Gk{uj) = 



(3.46) 



is the momentum-dependent Green's function. 

In general, the bare response function in Eq. (3.45) is 
a function of the frequency shift = cjj — cjj, of the 
incoming photon frequency uJi and the outgoing photon 
ojf frequency; it can be enhanced when one or both of 
the denominators are resonant (i.e., they coincide). In 
the latter case, we have a so-called "double" or "multiple 
resonance" .^^ The full response function also includes the 
vertex renormalizations. But the total (reducible) charge 
vertex in Eq. (3.21) for the Falicov-Kimball model does 
not diverge, and hence it does not introduce any addi- 
tional "resonances." It only leads to a renormalization of 
the total Raman response. 



IV. NUMERICAL RESULTS 

We begin our results by showing the single particle 
density of states of the spinless Falicov-Kimball model in 
infinite dimensions with (pe) = {tVi) = 1/2. The density 
of states is independent of temperature, and a metal- 
insulator transition occurs at f/ = \/2. In the insulating 
phase, the self energy develops a pole at w = 0, and the 
Green's function vanishes there. There is no true gap to 
this system, as the bare Gaussian density of states forces 
the interacting density of states to be nonzero whenever 
the self energy is finite. In Fig. 4, we plot the DOS 
for 5 values of U ranging from a weakly scattering metal 
U = 0.5, to a strongly scattering metal U = 1, to a near- 
critical insulator U = 1.5, a "small-gap" -insulator U = 2 
and a "moderate-gap" -insulator [7 = 3. Note that the 
metal-insulator transition is continuous for the Falicov- 
Kimball model, in the sense that the zero-temperature dc 
conductivity continuously goes to zero at the transition. 
Note further that in the metallic phase, the system is not 
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Frequency (w/t*) 

FIG. 4: Interacting singlc-particlc density of states for 
U = 0.5, 1.0, 1.5, 2.0, and 3.0 {U increases as the pseudo- 
gap becomes stronger). Note how the DOS first develops a 
depression near the chemical potential and then develops a 
pseudogap as the metal-insulator transition occurs (the DOS 
vanishes only at w = in the insulator). 



a Fermi liquid because the scattering time at the putative 
Fermi surface does not become infinite as T ^ 0. 

Once the self energy and the DOS are known, the dif- 
ferent contributions to Raman scattering can be deter- 
mined by straightforward, but tedious numerical integra- 
tions of the relevant functions for each scattering channel 
[Eqs. (3.13, 3.14, 3.25, 3.26, and 3.36-3.44). There are 
some subtleties with this approach, especially in the in- 
sulating phase, as the iterative approach to determining 
the DOS and the self energy becomes inaccurate once the 
imaginary part of the self energy becomes smaller than 
about 10~^^. Fortunately, there is a simple analytic form 
that can be used to construct the imaginary parts of the 
Green's functions and self energies in this regime, so all 
relevant quantities can be evaluated with care.^° 

We find that the Stokes response is significantly larger 
than the anti-Stokes response in the resonant regime, be- 
cause the double resonance greatly enhances the signal 
when the transfered energy approaches the incident pho- 
ton frequency (in the nonresonant regime, both Stokes 
and anti-Stokes responses are identical). Hence, we will 
present only the Stokes response here. We also find that, 
generically, the response "sharpens" as T — > 0, with the 
spectral response growing at low temperature (except for 
the low-energy, thermally excited response in the insulat- 
ing phase). Thus, we focus on low and moderate temper- 
atures in the metallic regime, since the Raman signal is 
largest there. 

In Figs. 5 and 6 we plot the total Raman response for 
U = 0.5 and U = 1 respectively. The former case is of a 
dirty metal, while the latter case is a metal that has such 

strong scattering that the density of states is depressed 
near the Fermi energy (but not so much as to create an 
insulator). The Stokes branch of the Raman response 
behaves in many respects as expected. The double res- 
onance causes a large enhancement of the signal as the 
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FIG. 5: Stokes Raman response for the three symmetry chan- 
nels in a dirty metal with U = 0.5. The Raman scattering 
response function is plotted as a function of the transfered 
frequency for incident photon frequencies ranging from 0.25 
to 4.5 in steps of 0.25 (the thickness of the lines aids in distin- 
guishing the different curves) . This data is at low temperature 
(T = 0.05) where the results are the "sharpest" . 



FIG. 6: Stokes Raman response for the three symmetry chan- 
nels in a strongly scattering "metal" with U = 1.0. The Ra- 
man scattering response function is plotted as a function of 
the transfered frequency for incident photon frequencies rang- 
ing from 0.25 to 4.5 in steps of 0.25. This data is at a moderate 
temperature (T = 0.5) where the nonresonant response has 
enhanced low-energy spectral weight in the Big channel. 



transfered frequency approaches the incident photon fre- 
quency. In the Loudon-Fleury regime, where the photon 
energy is much larger than the band energies, one can see 
a nice separation of the signal into the nonresonant and 
resonant (plus mixed) pieces (note that the nonresonant 
Aig response is small due to screening effects and the non- 
resonant i?2g response vanishes due to symmetry, but the 
resonant effects are strong in both of these channels). In 
general, the resonant effects are strongest near the dou- 
ble resonance, and it is not true that the total response 
looks like the nonresonant response plus a uniform res- 
onant enhancement, so resonant effects must be studied 
with care to understand the effects they play on the light 
scattering. Finally, note the overall similarity between 
panels (b) and (c) in Figs. 5 and 6. This arises from the 
fact that generically, the resonant effects overwhelm both 
nonresonant effects and mixed scattering effects, and it 
shows that there is not a huge variation in the resonant 
Raman response due to the additional renormalizations 
in the Aig channel. 

The insulating phase {U > \/2) provides a number of 
interesting new features to the electronic Raman scatter- 
ing (results for the near-critical insulator^^ at U = 1.5 
and for the small-gap insulator^*^ at J7 = 2 have already 
appeared). We begin with a discussion of a good cor- 
related insulator U = 3, which appears to have a well 
defined gap region in the DOS (but note that the DOS 
only vanishes exactly at w = 0). Hence we expect there to 
be significant thermally driven effects in this case. To be- 
gin, we plot the Raman scattering at a fixed temperature, 
but with varying incident photon frequency in Fig. 7. 



Note that there is substantial spectral weight in both a 
low-energy and a high-energy peak, and that when the 
incident photon frequency is approximately equal to U, 
the high-energy (charge-transfer) peak can be enhanced 
significantly. But something strange occurs for higher 
frequencies in the Aig channel. As increases beyond 
about 3.25, we stop to see the development of a separate 
charge transfer peak, and the net scattering curve looks 
like a simple double resonance curve even though the 
nonresonant response has a well-developed charge trans- 
fer peak. In other words, we are not seeing the evolution 
of the scattering to a simple break up of a nonresonant 
piece and a double resonance piece as Wi is made large. 
This may not be too surprising, because in the Aig chan- 
nel we have nonresonant, resonant, and mixed contribu- 
tions to the scattering. To illustrate how this occurs, we 
plot the separate contributions to the Raman scattering 
in Fig. 8 for the Big and Aig channels for Wj = 4.0. In the 
top panel, we see the expected shape for the nonresonant 
curve, with both low and high energy peaks, but sur- 
prisingly, there is a strong resonant enhancement of both 
peaks. This is even more dramatic in the bottom panel, 
where the vertex corrections suppress the nonresonant 
low-energy peak in the Aig channel, but the resonant 
terms bring back a strong enhancement in that region 
(in essence because the conservation of total charge acts 
to effectively screen the low-energy excitations, but the 
screening is much less effective for the resonant terms). 
The mixed contribution is small at low energy, but has 
a well developed charge-transfer-like feature, that is neg- 
ative, and completely overwhelms, and cancels the non- 
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FIG. 7: Stokes Raman response for the three synnnetry chan- 
nels in a correlated insulator with U = 3.0. The Raman scat- 
tering response function is plotted as a function of the trans- 
fered frequency for incident photon frequencies ranging from 
0.25 to 5.0 in steps of 0.25. This data is at a high tempera- 
ture (T = 1.0) where the nonresonant response has enhanced 
low-energy spectral weight in Big and Aig channels. 



resonant charge-transfer peak, leaving behind essentially 
a double rcsonancc-likc curve. These results are obvi- 
ously quite complex. If the incident photon frequency 
increases further, then the peak in the mixed response 
moves to higher energy, and the nonresonant peak plus 
a higher frequency double resonance peak picture holds, 
but the width of the double resonance peak can be ex- 
tremely narrow. One might be surprised that the double 
resonance peak survives in the insulator (because there 
are no electronic states within the gap), but in this case, 
we only have a pseudogap, and the states in the "gap 
region" are few in number, but long-lived and hence con- 
tribute to the scattering.^^ 

One of the common features in resonant Raman scat- 
tering is a large enhancement of the scattering when a 
new scattering channel opens, as the photon frequency 
becomes larger than an energy gap, for example. One 
question to ask is does such a feature survive in a corre- 
lated system. As described above, there is no energy gap 
in the insulating phase (on the hypercubic lattice), but 
there is a region where the DOS is exponentially small, 
and then increases rapidly to be of order unity. One can 
ask whether there are features in the Raman scattering 
that show enhancements when the photon frequency is 
larger than the width of the exponentially small "gap re- 
gion" of the DOS. Since the gap region for U = 3 is about 
0.5 above and below the chemical potential, we expect 
interesting results for photon frequencies near 0.5. We 
plot the Raman scattering for iVi increasing from 0.1 to 
1 in steps of 0.1 in Fig. 9 for low temperature (T = 0.2). 
Note how small the overall scale of the Raman scatter- 
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FIG. 8: Separation of different contributions to the Stokes re- 
sponse for U — 3,T = 1.0, and cui — 4.0 in the (a) Big and (b) 
Aig channels. The solid line is the total response, the dotted 
line is the resonant piece, the dashed line is the nonresonant 
piece, and the chain-dotted line is the mixed contribution. 



ing is. We see different behavior in the Big and i?2g 
sectors versus the Aig sector. In panels (a) and (b) we 
see the low energy scattering increases as uji increases 
until reaches approximately 0.5, where it starts to de- 
crease. The increasing behavior is essentially this reso- 
nant enhancement due to the opening of scattering chan- 
nels as the photon frequency becomes larger than the 
gap. Note how this phenomenon essentially does not oc- 
cur in panel (c), where the curves lie below each other 
as cJi is increased. Hence the Aig channel does not show 
the analogue of this resonant-enhancement effect. The 
effect disappears in all channels once the temperature 
becomes larger than about 0.5, where thermal excita- 
tions can be easily made across the "gap region". Note, 
furthermore, that the largest resonant effects occur not 
when the scattering channel first opens, but rather when 
LUi U because that is the value of frequency that sepa- 
rates the peaks in the single-particle DOS, and hence it 
corresponds to the strongest scattering from occupied to 
unoccupied states. 

We saw in Fig. 8 that there is a resonant enhancement 
at low energy when the incident photon frequency is close 
to U in size. To examine this phenomenon further, we 
plot the total Raman scattering at a fixed transfered pho- 
ton frequency (chosen to be 0.5 for the low-energy peak 
and 3.0 for the high-energy peak) as a function of the 
incident photon frequency in Figs. 10 and 11. 

In Fig. 10, we see expected behavior. The charge- 
transfer peak at n = [/ = 3 has a resonant enhancement 
for photon frequencies slightly higher than U, then a sup- 
pression to the nonresonant peak values at the highest 
incident frequencies (except for the Aig channel, where 
the charge transfer peak is initially suppressed until the 
incident photon frequency is larger than about 6, due 
to the cancellation from the mixed diagrams described 
above). The width of the resonant peak is about 0.5, 
and it is pushed to higher frequency in the Big and B2g 
channels. In Fig. 11, we find an interesting joint res- 
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FIG. 9: Raman response at low energy for (7 = 3 and T — 0.2. 
The incident photon frequency changes from 0.1 to 1.0 in steps 
of 0.1. 
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FIG. 10: Raman response at f2 = 3 for J7 = 3 and various 
temperatures. The horizontal axis is the incident photon fre- 
quency. The thickest curve is T = 0.05, and the temperature 
increases to 0.2, 0.5, and 1 as the curves are made thinner. 
Note that the curves for the lowest two temperatures are the 
largest, and are hard to separate. 



onancc effect. There is a resonant enhancement near 
uji = 0.5, that comes from the double resonance. In ad- 
dition, there is another broad resonance effect centered 
just slightly higher than Ui = U = 3, where both the 
charge-transfer and the low-energy peaks resonate at the 
same incident photon frequency. In the Aig channel, the 
joint resonance peak is a single smooth peak, while in the 
Big and channels, the joint resonance peak seems to 
have a double-peak structure to it. As the temperature 
is reduced, the resonant effects remain, but the spectral 
weight in the low-energy peak gets suppressed to very 
small values (the T = 0.05 curves are indistinguishable 



FIG. 11: Raman response at f2 = 0.5 for U — 3 and various 
temperatures. The horizontal axis is the incident photon fre- 
quency. The thickest curve is T = 0.05, and the temperature 
increases to 0.2, 0.5, and 1 as the curves are made thinner. 
Note that the T = 0.05 curves are more than three orders of 
magnitude smaller than the T = 0.2 curves, and cannot be 
distinguished from the horizontal axis in the figure. 



from the horizontal axis because they are at least three 
orders of magnitude smaller than the T ~ 0.2 curves). 
The evolution of the resonant profile for other values of 
transferred frequency O is complex and can be found in 
Ref. 25 for [7 = 3. 
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FIG. 12: Low-energy isosbestic behavior of the resonant 
Raman response for f/ = 3 and uoi — 3.5. The different 
thicknesses correspond to different temperatures, with thicker 
curves corresponding to lower temperatures (four curves are 
plotted for T = 1, 0.5, 0.2, and 0.05). Note how the curves all 
cross at an isosbestic point, just slightly smaller than U/2. 



The low-energy isosbestic behavior (which means that 
the Raman response is independent of temperature at a 
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characteristic frequency) is plotted in Fig. 12. We choose 
ijji = 3.5 because it corresponds to the maximal joint res- 
onance effect for both the charge-transfer and low-energy 
peaks. We find that the low-energy isosbestic behavior 
is generic for the resonant Raman scattering, with the 
response curves crossing at O k, U/2 = 1.5 for all sym- 
metries. Hence, the low-energy isosbestic behavior seen 
in the nonresonant response (which was most apparent in 
the i?ig channel, but can also be seen in the Aig channel 
when the response is plotted on a logarithmic scale^^), 
survives in the resonant cases as well, and this helps ex- 
plain why it is seen in so many experimental systems. 
In addition to the low-energy isosbestic point shown in 
Fig. 12 at f2 « U/2, there is a second isosbestic point^^ 
that appear near the double resonance O « Wj. Start- 
ing from large u)i (fii > U), as is reduced, the two 
isosbestic points move closer to each other, eventually 
joining together and disappearing when w C//2. So 
the isosbestic behavior will not be seen if the incident 
photon frequency is too low. 



V. DISCUSSION 

With the use of DMFT, we solved for the full Ra- 
man response for all frequencies of incoming light in the 
Falicov-Kimball model. Since the Falicov-Kimball model 
can be tuned across a metal-insulator transition, we have 
determined the form of Raman scattering in both the 
metallic and insulating states, and have investigated light 
scattering on both sides of the quantum critical point at 
U = \/2. Resonant, non- resonant, and mixed contribu- 
tions have all been treated on an equal footing and we 
allowed for an analysis of the dependence of Raman scat- 
tering with temperature, interactions, and different light 
polarizations. 

Our results confirm a number of previously held be- 
liefs. First, we find a strong resonant enhancement of the 
charge-transfer peak in Raman scattering when the inci- 
dent photon energy lies near the charge-transfer energy. 
This behavior is robust to temperature and polarization 
changes due to the local nature of the charge-transfer ex- 
citation in our model. Second, we also find a polarization- 
independent "double-resonance" enhancement when the 
transfered frequency of the light approaches the incident 
light frequency. This feature survives in the insulating 
phase because of the pseudogap nature of the insulator 
on the hypercubic lattice. 

In addition, we find a number of new features of light 
scattering in correlated insulators. We find that low en- 
ergy spectral features, related to thermal populations of 
elementary excitations, show resonance behavior when 
the incident light is tuned to the much higher frequency 
of the charge-transfer energy. This is a specific case where 
the correlations are crucial, since in uncorrelated mate- 
rials, this would correspond to off-resonant conditions. 
Yet due to the many-body nature of the correlated band, 
spectral features well separated from the charge trans- 



fer peak have a non-trivial resonance profile. We believe 
that these may be potentially useful to understand the 
complex nature of charge excitations in correlated ma- 
terials as it would impact both electronic and phononic 
Raman scattering at low frequencies. Finally, we find 
that the presence of an isosbestic point in the Raman re- 
sponse for correlated insulators results from a symmetry- 
dependent combination of all resonant, mixed, and non- 
resonant terms, and appears to be generic. 

We close with a discussion of open questions con- 
cerning improvements to the theory. Here we have re- 
stricted ourselves to Raman scattering in a correlated 
band of electrons in the limit of large spatial dimensions. 
Performing calculations in physical dimensions requires 
more many-particle charge vertex renormalizations which 
makes the problem extremely difficult, though possible 
in principle. But we found that most vertex renormal- 
izations were rather mild, so including nonlocal effects 
into the vertices (finite dimensions) probably does not 
change these results dramatically (unless the vertex can 
diverge in finite dimensions). In addition, q-dependent 
information would prove to be useful for investigating dis- 
persive many-particle excitations, as probed in inelastic 
x-ray scattering. These are topics of future interest. 
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APPENDIX A: PARQUET CONTRIBUTIONS 

In addition to the diagrams presented in Fig 3, there 
can also be parquet-like contributions with both vertical 
and horizontal renormalizations. One type of these dia- 
grams is shown in Fig. 13. The corresponding expression 

has the form 

qkk' 

al3a'l)' 

(Al) 

xr2r(l,3)r(2,4)Gi(fc')G2(fc')G2(fe' + q)G^{k' + q) 
xG2{k + g)G4(fc + q)G^{k)Gi{k). 

In the expression in Eq. (Al), all the momentum depen- 
dence (in the D = oo limit) is contained in the band 

energy [see Eq. (3.3)] and, after expanding the products 
of the Green's functions with the same momentum into 
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partial fractions over Cfe, the summations over momen- 
tum are of the form 



TV 



E 



1 sin fca sin(fc^ + g^) 



(A2) 



sinfc^, sin(fc^, + 
^j^t^, (^3 - fife') (^4 - ek'+q) 
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k+qi ^^ \^\ J k'+q 



In the same way as was done for Eq. (3.16), we find 
the expression in the bracket reduces to the following in 
the Z) — > 00 limit [Q = sgn(ImZi)] 

Cioo C200 


(A3) 



i^COSga 



AD 



sin sin 



where Xq = lim£)_+oc X]a=i '^'^^'Za- The main contri- 
bution comes from the first term in Eq. (A3) and there 
is a similar term in the second bracket of Eq. (A2). Re- 
placing the square of a cosine by its average value 5 , we 
find that Eq. (A2) reduces to 



lim GiG2C?3G4^0 

D->cx) 8D 



(A4) 



FIG. 13: Feynman diagrams for a typical paxquet-like renor- 
malization. This resonant diagram has a simultaneous hor- 
izontal and vertical renormalization by the two-particle re- 
ducible charge vertex. Note that such a renormalization is 
only possible in the sector. 



which vanishes as ^ 00. A similar procedure can be 
performed for all other terms with a parquet-like renor- 
malization. Hence, the parquet-like contributions are 
unimportant for resonant Raman scattering in large di- 
mensions. 
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